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Abstract 

Chloroplasts originated from cyanobacteria through endosymbiosis. The original cyanobacterial endosymbiont evolved to 
adapt to the biochemically rich intracellular environment of the host cell while maintaining its photosynthetic function; 
however, no such process has been experimentally demonstrated. Here, we show the adaptation of a model 
cyanobacterium, Synechocystis sp. PCC 6803, to a biochemically rich environment by experimental evolution. Synechocystis 
sp. PCC 6803 does not grow in a biochemically rich, chemically defined medium because several amino acids are toxic to 
the cells at approximately 1 mM. We cultured the cyanobacteria in media with the toxic amino acids at 0.1 mM, then serially 
transferred the culture, gradually increasing the concentration of the toxic amino acids. The cells evolved to show 
approximately the same specific growth rate in media with 0 and 1 mM of the toxic amino acid in approximately 84 
generations and evolved to grow faster in the media with 1 mM than in the media with 0 mM in approximately 181 
generations. We did not detect a statistically significant decrease in the autotrophic growth of the evolved strain in an 
inorganic medium, indicating the maintenance of the photosynthetic function. Whole-genome resequencing revealed 
changes in the genes related to the cell membrane and the carboxysome. Moreover, we quantitatively analyzed the 
evolutionary changes by using simple mathematical models, which evaluated the evolution as an increase in the half- 
maximal inhibitory concentration (IC 50 ) and estimated quantitative characteristics of the evolutionary process. Our results 
clearly demonstrate not only the potential of a model cyanobacterium to adapt to a biochemically rich environment 
without a significant decrease in photosynthetic function but also the properties of its evolutionary process, which sheds 
light of the evolution of chloroplasts at the initial stage. 
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Introduction 

Chloroplasts originated from cyanobacteria through endosym- 
biosis [1-4]. Chloroplasts conduct photosynthesis and synthesize 
most fatty acids and some amino acids and hormones [5]. 
Although chloroplasts play many important roles in plant cells, 
they lack many genes and functions compared to cyanobacteria, 
depending instead on the host plant cells. These data suggest that 
the original cyanobacterial endosymbiont adapted to the environ- 
ment in plant cells, which is biochemically richer than the 
environment of free-living cyanobacteria [6]. Thus, the original 
cyanobacterial endosymbiont evolved to depend on biochemical 
substances provided by the host plant cell, maintaining their 



photosynthetic functions to provide that cell with biochemical 
substances in return. 

It is not clear how cyanobacteria can adapt to a biochemically 
rich environment while maintaining their photosynthetic function. 
Cyanobacteria inhabit almost every niche on earth, including 
freshwater and salty, hot, and arid environments [6], and employ 
various stress-management mechanisms [7]. Cyanobacteria also 
establish long-lived symbiotic associations with fungi and within 
the plant kingdom [8]; however, cyanobacteria are sensitive to 
many important biochemical substances. For example, in the case 
of Synechocystis sp. PCC 6803, a prominent model cyanobacteria, 
the original strain (stocked in the Pasteur Culture Collection) is 
sensitive to glucose in light [9]. While glucose-tolerant lines exist, 
they are still sensitive to several amino acids [10-12], possibly 
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because their autotrophic ability opposes the ability to use nutrient 
biochemical substances. Indeed, Bell and Reboud [13-15] 
experimentally evolved Chlamydomonas reinhardtii in light (photoau- 
totrophic) and dark (heterotrophic) environments and showed that 
some lines in the dark environment evolved to have higher specific 
growth rates in the dark environment but lower specific growth 
rates in the light environment, although C. reinhardtii is not a 
cyanobacterium. The experimental demonstration showing that 
cyanobacteria adapt to a biochemically rich environment sheds 
light on the mechanism of adaptation and how cyanobacteria 
evolved to become chloroplasts. 

In this study, we experimentally evolved Synechocystis sp. PCC 
6803 to adapt to a biochemically rich environment without loss of 
autotrophic ability. We used a synthetic medium (TCM1) that is 
sufficiently nutrient-rich to allow the growth of a heterotrophic 
model protozoan, Tetrahymena thermophila, and to inhibit growth of 
the cyanobacteria. We first identified the six amino acids in TCM 1 
that are toxic to the cyanobacteria (Figure SI in File SI). We then 
cultivated the cyanobacteria in the medium with reduced 
concentrations of the toxic amino acids and serially transferred 
the culture, gradually increasing the amino acids concentrations. 
We found that the cyanobacterium adapted to the rich environ- 
ment in approximately 100 generations, and we detected no 
significant difference in autotrophic growth in an inorganic 
medium between the ancestral and evolved populations. Whole- 
genome resequencing revealed changes in genes related to the cell 
membrane and to the carboxysome, which is involved in 
photosynthesis. Our results clearly show that the model cyano- 
bacterium can evolve to adapt to a rich environment with no 
significant decrease in photosynthetic function and experimentally 
show the properties of the evolutionary process of the cyanobac- 
teria to adapt to a biochemically rich environment. 

Materials and Methods 

Strains and culture conditions 

The Synechocystis sp. PCC 6803 (GT) strain, originally obtained 
from Dr. J. G. K. Williams (Dupont de Nemours, Wilmington, 
DE, U.S.A.), was used as the ancestral strain. This is a glucose- 
tolerant strain and not the original strain (stocked in the Pasteur 
Culture Collection). All cells were grown at 30°C in 10 mL liquid 
culture in a 50 mL Erlenmeyer flask rotated at 80 rpm using NR-3 
(Taitec, Saitama, Japan), with 10 |jmol-m -2 -s _1 light using an 
LED system (MIL-C1000T, MIL-U200, MIL-R18A, and MIL- 
B18A; Sanyo, Osaka, Japan). The LED board ratio of red (MIL- 
R18A) to blue (MIL-B18A) was 9. We used TCM1, TCM0, BG- 
11, and BG-ll+AA for the culture media. TCM1 was derived 
from a chemically defined medium for a heterotrophic model 
protozoan, Tetrahymena thermophila CDM15 [16], with a few 
modifications (see Table SI in File SI for the detailed compo- 
nents). TCM1 includes the following 6 amino acids that are toxic 
to cyanobacteria: Arg, His, Lys, Met, Phe, and Thr in 
concentrations of approximately 1 mM. TCM0 does not include 
the 6 toxic amino acids and the other components are the same as 
the components in TCM1. We also used TCMx, where x is the 
concentration of the 6 toxic amino acids relative to TCM1 (e.g., 
TCM0.1 when x = 0.l). BG-11 is an inorganic medium that is 
standard for Synechocystis sp. PCC 6803 culture [17]. BG-ll+AA 
includes the 6 toxic amino acids in the same concentrations as the 
TCM 1 and the other components are the same as BG-11. 

Evolution Experiment 

We first cultivated the cells (pre-cultured in BG-1 1) in TCM0, 
transferred them 4 times (as an initial adaptation, see Figure S2 in 



File SI), freeze-stocked them once by adding DMSO to the culture 
(final 7%), and recovered them in TCM0 before beginning the 
experimental evolution. The culture conditions were the same as 
in the evolution experiments except for the culture media. Simply, 
the difference between BG-1 1 and TCM0 is that TCM0 includes 
amino acids, vitamin, nucleosides, and glucose but BG-11 does 
not. We then serially transferred the cell culture following the 
protocol described below, always under 10 umol-m _2 -s _1 of light. 
At every transfer until day 79, we split the cell culture into three or 
four TCMjc media that differed in x. After the cultivation, we 
identified the culture with the highest x that showed a population 
increase of more than the square root of 1 0, and we used it as the 
origin for the next transfer. The initial cell concentration of each 
transfer was set to be approximately 2xl0 6 cells/mL (thus the 
dilution rate for transfers was different for each round and the 
bottleneck size of the population was 2x10 cells in 1 0 mL), and 
the actual cell concentrations in the evolution experiment are 
plotted in Figure 1A. The incubation time for each culture was 3 
or 4 days so that we maintained almost exponential phase for the 
transferred culture (under approximately 2x10 cells/ mL) as semi- 
continuous culture. After day 79, we serially transferred the cell 
culture to TCM 1 , examining the growth in TCM0 in parallel until 
day 271 (see Figure S4 in File SI). The generation during 
evolution was calculated as ln(A rl , op )/ln2, where is the total 

fold change in the population. 

Measurement of cell concentrations 

We measured the cell concentrations relative to a known 
concentration of fluorescent beads (Fluoresbrite YG Microspheres, 
6 um; Polysciences Inc., Warrington, PA, USA) using a Cytomics 
TM FC500 flow cytometer (Beckman Coulter, Inc., CA, USA) by 




Culture time (day) Culture time (day) 

Figure 1. Evolutionary changes during adaptation to a 
nutrient-rich environment. We cultivated the cells in media with 
various concentrations of toxic amino acids until culture day 79, then 
used only TCM1 and TCM0. (A) Population dynamics in the evolution 
experiment until day 79. The colors show the culture media used (TCMx 
media; x is shown on the right). The transfers of the culture by dilution 
are depicted as the vertical decrease in the cell concentration. (B) 
Evolutionary changes in the specific growth rate (/<) in media with 
different toxic amino acid concentrations (the colors are same as in A). 
The specific growth rate was determined by the slope of the linear 
regression of the natural log of cell concentration per day, using first 
three (or two, if only two points were measured in the culture) data 
points. (C) The evolutionary changes in the relative specific growth rates 
of the cells in TCM1 (/(,) and in TCM0 (/i 0 ). The raw data of the 
population dynamics are shown in Figure S4 in File S1. 
doi:10.1371/joumal.pone.0098337.g001 
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loading culture samples mixed with the beads as described 
previously [18]. A 488-nm argon excitation laser was employed, 
and bandpass niters of 515-535 and 610-630 nm were used to 
measure green and red fluorescence, respectively. The cyanobac- 
terial cells were detected as red fluorescent particles due to the 
presence of chlorophyll, and the cell concentration was calculated 
from the counts of these particles relative to the fluorescent beads, 
which were detected as green fluorescent particles. 

Whole-genome resequencing 

We analyzed 3 genome sequences of (i) ancestral cells, (ii) 
initially adaptive cells in TCMO (Figure S2 in File SI), and (iii) 
evolved cells (at 408 days, 85 transfers, 548 generations, in 
TCM1). We prepared precultures by shaking stocked populations 
in media of BG-11, TCMO, and TCM1 for (i), (ii), and (iii), 
respectively, under the same culture conditions as described above, 
until the cell concentrations reached approximately 1 xlO 8 cells/ 
mL. The cyanobacterial cells were collected by centrifugation at 
14, 170 xg for 2 min at 25°C, and the pelleted cells were stored at 
— 80°C prior to genomic DNA purification. Genomic DNA was 
isolated and purified in accordance with a standard method [19]. 
Genome sequence analyses were performed with an Illumina 
MiSeq sequencer (Illumina, CA, USA) by Eurofms MWG Operon 
(AL, USA). Three hundred basepair insert libraries were 
generated according the Illumina protocol and sequenced in an 
Illumina MiSeq. In this study, 3 samples with different indexes 
were mixed and then sequenced in 2x150 bp sequencing mode. 
We obtained 1079, 917, and 1110 Mbp sequences for (i), (ii), and 
(iii), respectively. These numbers of base calls corresponded to 
x301.9, x256.6, and x310.6 redundancy, respectively. We used 
Iuurnina's CASAVA pipeline version 1.8.2 software (Illumina, CA, 
USA) for data analyses including trimming and alignment to the 
reference genome sequence NC_000911.1 (Synechocystis sp. PCC 
6803 chromosome, complete genome). We used VarScan 2 
software [20] for variant analysis. 

Modeling and simulations 

Simulations of the evolutionary change in lnflCso] were 
performed using Matlab software (MathWorks, MA, USA). The 
Matlab code used is presented in Appendix S 1 in File S 1 . 

Results and Discussion 

Evolutionary changes in adaptation to a nutrient-rich 
environment 

We serially transferred the cyanobacterial cell culture, gradually 
increasing the concentrations of the 6 amino acids that were toxic 
to the cells and always maintaining the cultures under light to 
maintain their photosynthetic function. Figure 1A shows cell 
growth during evolution; the cells evolved the ability to grow in 
media with higher concentrations of the toxic amino acids through 
serial transfer. Figure IB shows the specific growth rates (jl) in 
media with different toxic amino acid concentrations, TCMa-, 
where x is the relative amino acid concentration (plots of fi as a 
function of x are shown in Figure S3 in File SI). We found a 
statistically significant increase in ji in TCM1 red line) while 
we did not find statistically significant increases in \l in TCMO (/lo, 
black line) by a Student's /-test to compare the slope of the linear 
regression of \i on the culture time with zero slope using a 
significance level of 0.05; the slopes were 0.0069 (7^ = 0.32, the 
degree of freedom DF= 20, p = 0.0057) and -0.0014 (it 2 = 0.032, 
DF= 19, /> = 0.44)/day 2 for TCM1 and TCMO, respectively, p, in 
TCMa for 0<,t<l (colored lines) were in between fii and fa. 
These results suggest that the key change was related not to global 



growth but to tolerance of the toxic amino acids. Then, we 
continued the transfer in TCM1 (the population dynamics are 
shown in Figure S4 in File SI). Figure 1C shows the evolutionary 
changes in fi t relative to fi Q . [i\ approached 90% of fig at 
approximately 75 days (2 1 transfers, 84 generations) and exceeded 
jUo at approximately 160 days (40 transfers, 181 generations). 
These results indicate that the model cyanobacteria had adapted 
to a biologically rich environment, increasing their tolerance to the 
previously toxic amino acids. 

Autotrophic ability of the evolved strain 

We found no statistically significant decrease in autotrophic 
ability (photosynthetic function) in the evolved cells. Figure 2 
shows the fi of the ancestral and evolved cells in various culture 
media. First of all, we found that the evolved cells actually 
increased their growth ability in TCM1 because there was a 
statistically significant difference in the mean fi in TCM1 between 
the evolved and the ancestral cells as determined by a Welch's t- 
test using a significance level of 0.05 (DF=7, p = 1.9xl0~ ). 
Importantly, there were no statistically significant differences 
between the evolved and the ancestral cells in the inorganic 
medium BG-11 (DF=5, p = 0.67), suggesting that the photosyn- 
thetic function was maintained in the evolved cells. Furthermore, 
the evolved cells showed amino acid tolerance even in BG- 1 1 
(there were no statistically significant differences between BG- 1 1 
and BG-ll+AA; DF=2, p = 0.U), while the ancestral strain was 
sensitive (there was a statistically significant difference between 
BG-11 and BG-ll+AA; DF= 4, p = 5.9 x 10" 6 ). We did not 
observe decreased means of the specific growth rates in the 
evolved cells (trade-offs) in any of the tested media. 

Whole-genome resequencing 

The evolved changes in phenotypes were corroborated by the 
genetic changes detected by whole-genome resequencing. We 
analyzed 3 genomic DNA samples of ancestral cells, initially 
adaptive cells in TCMO (Figure S2 in File SI), and evolved cells in 
TCM1. All mutations in the chromosomes found in the variant 
analysis are listed in Table S2 in File SI. We found 20 single-base 
substitutions and 8 small indels in ancestral cells compared with 
the reference sequence (NC_000911.1 Synechocystis sp. PCC 6803 
chromosome, complete genome). These mutations were also 
detected in the other two types (initially adaptive and evolved) of 
cells. The detected mutations were identical for the ancestral and 
initially adaptive cells (see Table S2 in File SI). 

The mutations detected only in the evolved strains are 
summarized in Table 1. All of the single-base substitutions 
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Figure 2. Specific growth rates of the ancestral and evolved 
cells in various environments. The black and white bars indicate the 
specific growth rates of the ancestral and evolved (at 408 days, 85 
transfers, 548 generations) cells, respectively, in the culture media 
shown at the bottom. The error bars show the SD of five (for BG-1 1 and 
TCM1 ) or two (for the other media) independent cultures. The ancestral 
culture in BG11+AA showed almost zero growth. 
doi:10.1371/joumal.pone.0098337.g002 
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Table 1. Genetic changes detected in the evolved cells relative to the ancestral cells. 





# 


Position (nt) 


Nucleotide 
changes 


Gene 


AA changes 
in the Gene 


1 


335496 


C to A 


Carboxysome-formation protein CcmA; complementary 


G to E 


2 


829508 


C toT 


RNA polymerase alpha subunit; complementary 


R to Q 


3 


1114921 


G to C 


Periplasmic substrate-binding and integral membrane protein of the ABC-type 
Bgt permease for basic amino acids and glutamine BgtB; complementary 


F to L 


4 


1128135 


C to G 


Unknown protein (sill 265); complementary 


A to P 


5 


Next to 1905171 


Del. GCCTCG 


Penicillin-binding protein (cell-wall biogenesis); complementary 


Del. AE 


6 


3203715 


G to A 


Probable cation transporter; complementary 


P to S 



All detected mutations in the chromosome of the ancestral cells, initially adaptive cells, and the evolved cells, relative to the reference (NC_000911.1; Synechocystis sp. 
PCC 6803 chromosome, complete genome), are shown in Table S2 in File SI. 
doi:1 0.1 371/journal.pone.0098337.t001 



identified in annotated genes were non-synonymous, suggesting 
that they were adaptive. Mutations in genes related to membrane 
biogenesis and transportations were detected (#3, 5 and 6 in 
Table 1), suggesting that the amino acid tolerance arose from 
changes in membrane functions, such as increased membrane 
protection from penetration by amino acids. Although gene #4 
was annotated as an unknown protein, it has been reported that 
this gene is related to stress responses, such as butanol, hexane, salt 
stresses and nitrogen starvation [21]. These results are consistent 
with the increase in tolerance. A mutation in the carboxysome- 
formation protein CcmA [22] (#1) might indicate an optimization 
of photosynthetic activity in a biochemically rich environment 
because the C0 2 -concentration mechanism (CCM) of the 
carboxysome is important to meeting habitat requirements [7]. 
A mutation was also identified in the RNA polymerase alpha 
subunit (#2). Although it is difficult to estimate the direct 
involvement of this mutation in the evolved changes in the 
phenotypes, the mutation in the RNA polymerase alpha subunit 
could potentially provide some benefit because the initiation of 
transcription is the main determinant for the regulation of gene 
expression in bacteria. Note that those individual mutations would 
not be directiy correlated with acquired phenotype such as by the 
existence of epistasis [23,24] and epigenetic changes [25]. Further 
work, including construction of a mutant containing these 
mutations and analysis of phenotypes of these mutants, is needed 
to elucidate the relationship between genotype and phenotype. 
Also note that we analyzed mixed (not cloned) populations that 
were stocked directly from the evolution experiment, and some 
heterogeneous mutations were found (see Table S2 in File SI). 
However, it is possible that this heterogeneity did not arise from 
the heterogeneity of different cell groups because it is known that 
Synechocystis sp. PCC 6803 has a multicopy genome [26], and 
heterozygous mutations are commonly found in the genome of the 
Synechocystis sp. PCC 6803 [27]. Actually, the heterogeneic 
frequencies were maintained during the course of evolution, and 
the specific mutations in evolved cells shown in Table 1 were 
almost homogeneous (see Table S2 in File SI). 

Quantitative analyses of the evolutionary changes using 
simple mathematical models 

We analyzed the experimental results of the evolutionary 
changes using simple mathematical models to extract the 
quantitative properties of the evolution. Specifically, we estimated 
a key trait (IC 5U , see below) in the evolution, the trait variance in 
the population, and some properties of the fitness landscape 
regarding the key trait, such as the probability of the beneficial 



mutation per replication and the extent of the change in the key 
trait per beneficial mutation. We used the data through day 75 (77 
generations) where the growth rate of the cells in the transferred 
cultures was less than [/q. First, we estimated the maximum specific 
growth rate (jum.^) and the half-maximal inhibitory concentration 
(IC 5 n) of the cell population for every transfer round using the 
experimental data of the relationship between the specific growth 
rate and the amino acid concentration (Figure IB, and Figure S3 
in File SI). Specifically, we estimated those two parameters by 
fitting the experimental data to a basic noncompetitive-inhibition 
model /.t = J Um ax /(l+x/IC 50 ) [28] (the fitting results are shown in 
Figure S3 in File SI). Figure 3A-i summarizes the fitting results 
and clearly shows that the principal increase was in IC 30 , not ^im f „; 
there was a statistically significant increase in lnflCso] ( re d square) 
while not in ln[um ax ] (black circle) as determined by a Student's t- 
test to compare the slope of the linear regression of the parameters 
on generation with zero slope using a significance level of 0.05; the 
slope of lnpC 50 ] and InCumJ was 0.035 {R 2 = 0.67, DF= 20, 
/> = 3.6xl0~ 6 ) and -0.0014 (R 2 = 0.037, DF=20, /; = 0.39) per 
generation, respectively. Thus, it is possible to quantitatively 
understand this evolutionary change as the change in a trait, IC 50 , 
based on quantitative genetics as described below. 

We analyzed the evolutionary rate of the IC 50 , specifically the 
slope of ln[IC5o] (designated as Sj C50 ), to provide an interpretation 
of the experimental results, using a simple model of the directional 
selection based on quantitative genetics [29-32] . Let us assume a 
cell population in which each cell has a trait z ( = ln[ICr, 0 ]) and a 
specific growth rate fi(z), with the frequency J{z,t) at the total N{i.e., 
2_j{z,t) = N). When a cell with z reproduces an offspring through 
replication, the trait of the progeny becomes z+d z by a beneficial 
mutation with a constant probability r, or else keeps z (probability 
1— r), which is the same as the parent (Figure 3B), i.e., r is a 
beneficial mutation rate per generation and d e is an effect size of a 
single beneficial mutation on the trait z. Constant r and d z indicate 
that this model assumes a non-epistatic landscape [33]. The rate 
equation oij(z,t) is described as 

= (i -rMz)f{z4) + ni(z-d z )f{z-d-_j)-]i(t)AzJ) (1) 

where Jx{t) is the mean specific growth rate of the population (see 
Text SI in File SI for the derivation). The first and second terms 
on the right-hand side represent the increase in the population 
frequency of z by reproduction from the population of z without 
any beneficial mutation (a gray square fraction at z in Figure 3B) 
and by reproduction from the population of z~ d z with a beneficial 
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mutation (a red square fraction at z~ d z in Figure 3B), respectively. 
The third term represents the decrease in the population by a non- 
selective dilution to keep the total population size (JV) constant. The 
rate of the evolutionary change in the mean of z over generations 
in this model has the same meaning as iSicso if constant. To obtain 
a simple understanding of the value of Sicso, we used a set of 
rough approximations to obtain an approximate analytical 
solution of Sj C50 from the model (see Text SI in File SI). We 
found that the evolutionary rate became constant in this model 
and obtained its approximate analytical solution as 

S/C50*2(ln2)4= In (2) 

where a st is a standard deviation of the trait z (the width of the trait 
frequency distribution) when the evolutionary rate is constant. 
Equation 2 indicates that 5ic50 is a monotonically increasing 
function of d z , r, and jVfor the ranges d z >0, 0<r<l, and N>\. 
The difference between this approximate analytical solution 
Equation 2 and the numerical simulation according to Equation 
1 is shown in Figure S5 in File SI. Equation 2 means that the 
evolutionary rate is proportional to the variance (a sl ), known as 
Fisher's fundamental theorem of natural selection [34], and 
<r J( ~0.16 was obtained from the experimental value of Sicso 
( = 0.035) independendy of the model parameters (see Figures 3A- 
ii and 3A-iii for graphical presentations of the frequency 
distribution). 

Undetermined variables were only d z and r because jVwas given 
in the experimental condition. The dependency of >$ic50 on d z and 
r is shown in Figure 3C. Therefore, if one of the two variables is 
determined, the other can also be determined from the 
experimental results, and the evolution can be visually reconsti- 
tuted by the model, as shown by the blue line in Figure 3A-i and 
the corresponding population changes in Figure 3A-ii using (d B 
r) = (0.4, 0.00023) as an example. This value of d z = 0A (i.e., 5ic50 
increases 1.5-fold per beneficial mutation) is possible because the 
extent of the total increase in lnpC 30 ] that is required to satisfy 
fio~fi\ (A z ) seems to be approximately 2.5 (—2 to 0.5) based on 
Figure 3A-i, and the total number of beneficial mutations (n m ) can 
be 6 based on the number of detected mutations, resulting in 
d z = A z / n,„~0.4. Then, r should be approximately 0.0002 to 
explain the experimental results by this model with d z = 0.4 as 
above. With r~0.0002 (i.e., 0.02% of replications are beneficial), 
the fraction of beneficial mutation among all possible mutations 
can be estimated as 0.07 (i.e., 7% of mutations are beneficial) from 
the value of r divided by the mutation rate per genome per 
replication m^^, which is known to be approximately 0.003 in 
microbes [35]. Note that r is more sensitive than d z in those range 
(Equation 2 and Figure 3C), e.g., if the range of d z is 0.3—0.5, the 
range of r is approximately 10 6 ~10 3 . Also note that above 
analyses are rough estimation and validity is still unclear. For 
example, the total number of beneficial mutations n m can be less 
than 6, for reasons such as non-beneficial mutations in Table 1 or 
few mutations until day 75. In addition, n m also can be greater than 
6 because heritable trait changes can occur not only through 
mutations but also by other epigenetic factors [25]. Although the 
validity of the values and the model itself are still unclear in this 
study, the above information provides important insights into the 
evolutionary process of the cyanobacteria to adapt to a biochem- 
ically rich environment. 




Figure 3. Quantitative analyses of the evolutionary changes by 
simple mathematical models. (A) Evaluation of the evolutionary 
changes in the maximum specific growth rate (^m ax ) and the half- 
maximal inhibitory concentration (IC 50 ) until day 79. (i) Evolutionary 
changes in ln[um ax ] (black circle) and ln[IC 50 ] (red square), determined 
by the fitting shown in Fure S3 in File SI , with respect to the generation 
(g) of each round. The error bars show the standard error of the fitting 
estimates. The black and red lines show the linear regression of \n[fjm ax ] 
and ln[IC 50 ], respectively (-0.001 4g+0.28 and 0.035g-2.5, respectively). 
The blue bold line is the change in the mean of ln[IC 50 ] calculated by 
the numerical simulation of the mathematical model of Eq. 1, using (d z , 
r, N) = (0.4, 0.00023, 10 7 ), assuming a homogeneous population for the 
initial state. The source code of the simulation is shown in Appendix S1 
in File S1. (ii and iii) Simulated evolutionary changes in frequency 
distributions of the cell population calculated by the numerical 
simulation using (d z , r, N) of (0.4, 0.00023, 10 7 : the same as the blue 
line in A-i) and (0.1, 0.04, 10 7 ), respectively. The vertical axis is ln[IC 50 ], as 
in A-i. (B) A conceptual scheme of the mathematical model of Eq. 1. See 
the text for explanation. Simply, the population shifts toward greater 
ln[IC 50 ] with a beneficial mutation rate per generation (r) and an effect 
size of a single beneficial mutation on ln[lC 50 ] (d z ). The axes correspond 
to those of A-ii and A-iii. (C) The dependency of the evolutionary rate of 
ln[IC 50 ], i.e., the slope of the red line in A-i (5| C 5o), on the model 
parameters d z and r calculated from Eq. 2 using N= 10 7 . The gray bold 
line shows the experimentally determined value (S| C50 = 0.03 5/g). The 
black point shows the value used for the simulation shown in A-ii, 
which is discussed in the text. 
doi:10.1371/joumal.pone.0098337.g003 

Conclusions 

In this study, we show that a model cyanobacterium, Synechocystis 
sp. PCC 6803, can evolve to adapt to a biochemically rich 
environment without statistically significant losses of autotrophic 
ability. Our results thus show the potential of a cyanobacterium to 
evolve into a chloroplast-like endosymbiont. Furthermore, we 
showed some properties of the evolutionary changes such as an 
increase in tolerance to toxic amino acids rather than the global 
growth rate itself, the quantitative characteristics of the evolution- 
ary rate, which provides estimated variance of the tolerance in the 
population and clues regarding the evolutionary fitness landscape, 
and genomic analyses that identified the mutated genes and were 
consistent with the increase in tolerance. Additional generations 
would be needed to show more qualitatively distinct changes 
toward a chloroplast-like phenotype, such as loss of growth ability 
in inorganic media while maintaining photosynthetic function and 
optimization of photosynthesis for a biochemically rich environ- 
ment. Community-level function (shown in another cyanobacte- 
rium [36]) must be a key factor for loss of functions in the 
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evolution [37]. Therefore, ultimately, the cyanobacteria cells 
would have to be cultivated with or within another eukaryotic cell 
[38-40] for the complete experimental evolution toward chloro- 
plast-like endosymbiont. 

Supporting Information 

File SI Supporting information which contains the following. 
Figure SI: Toxic components in TCM1. Figure S2: Preculture for 
the evolution experiment. Figure S3: The specific growth rates as a 
function of the relative concentration of the toxic amino acids in 
each transfer round. Figure S4: Population dynamics in the 
evolution experiment. Figure S5: The difference between the 
approximate analytical solution (Equation 2) and the numerical 
simulation according to Equation 1. Table SI: Components of 
TGM1. Table S2: All detected mutations in the chromosome, 
relative to the reference (NC_0009 1 1 . 1 ; Synechocystis sp. PCC 
6803 chromosome, complete genome), by the genomic analysis. 
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